Evolution of a mass-less test scalar field on Boson Stars space-times 
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We numerically solve the massless test scalar field equation on the space-time background of bo- 
son stars and black holes. In order to do so, we use a numerical domain that contains future null 
infinity. We achieve this construction using a scri-fixing conformal compactification technique based 
on hyperboloidal constant mean curvature foliations of the space-time and solve the conformally 
invariant wave equation. We present two results: the scalar field shows oscillations of the quasinor- 
mal mode type found for black holes only for boson star configurations that are compact; and no 
signs of tail decay are found in the parameter space we explored. Even though our results do not 
correspond to the master equation of perturbations of boson star solutions, they indicate that the 
parameter space of boson stars as black hole mimickers is restricted to compact configurations. 
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I. INTRODUCTION 

Observations related to high energy events, indicate 
the existence of black hole candidates (BHCs). It is usu- 
ally assumed that black hole solutions are the models of 
BHCs. Nevertheless there are various nonvacuum solu- 
tions in general relativity that are being studied as alter- 
natives to black holes, sometimes called black hole mim- 
ickers. Among such alternative solutions are wormholes 
[l|, 0] , gravastars [l|, y-Q , branc world solutions Q and 
boson stars |7H10| . Predictions on phenomena related to 
black hole mimickers are important, either because they 
may rule out the mimicker or at least restrict the param- 
eter space of mimicker configurations. 

One of the most important properties of mimickers is 
the stability. For viformhole solutions instability has been 
shown for basic wormholc solutions |lll4l4| . destroying 
previous hope on the possibility that these are stable 
as shown in [15|. The stability of gravastars has been 
explored and it has been found that there are stability 
regions 0, [l6l - [l8t . In favor of boson stars the stabil- 
ity of the solutions has been exhaustively studied, using 
perturbative methods |19l . |20| and using full nonlinear 
numerical relativity, both in spherical symmetry j2ll |22| 
and full three dimensions [23|. 

In the particular case of boson stars they have been 
studied as sources of gravitational radiation, both as bi- 
nary systems Q and as perturbed boson stars [231 in the 
full nonlinear regime. The study of gravitational wave 
signatures has also been pursued in order to differentiate 
between a black hole and a gravastar [j]. Instead, in the 
case of wormholes for example, simple solutions (sup- 
ported by a phantom scalar field) would not stand the 
perturbative analysis and there is no hope for a study of 
a binary system, because as shown in [l3| the lifetime of 
these solutions is rather short and they should either col- 
lapse and form black holes or explode before they could 
merge (although in |14l | it was shown that charged worm- 
holes can have a longer lifetime if the charge parameter 
is adequately chosen). At the end of the day, if black 



hole mimickers are to be compared from all the angles 
with black holes, mimickers should also be expected to 
do what black holes can do: they collide and generate 
gravitational radiation with a given fingerprint. 

In this paper we explore a simple but potentially fruit- 
ful problem: the evolution of a test massless scalar field 
on the space-time of boson stars. It is well known that, in 
the case of black holes, the solution of such an equation 
is related to the quasinormal modes of axial perturba- 
tions and electromagnetic test fields [26[. And even if 
this equation is not the perturbation equation of boson 
stars (it would involve the perturbation of matter too) 
the behavior of the massless scalar field would provide 
indications on which boson stars look more like black 
holes, among the big set of boson stars that are mimickers 
when the power emitted by accretion disks are compared 
to those due to black holes [3, [l3| • We calculate the nu- 
merical solution for various boson star space-times and 
compare to the solution for the space-time corresponding 
to a Schwarzschild black hole. We track the numerical 
solution at the level of mode ringing and tail decay. 

One of the important ingredients of this paper is that 
we foliate our space-times with hyperboloidal slices, for 
which we develop the procedure to foliate a static spher- 
ically symmetric space-time with constant mean curva- 
ture slices. Following [23, [S^' "^^ compactify the bound- 
ary at j/+ (future null infinity) which allows one to mea- 
sure the amplitude of the scalar field at almost arbitrary 
distances from the source. The compactification of future 
null infinity implies the physical metric to be singular at 
infinity; this problem is fixed using a rescaling through 
a conformal transformation [27[. Therefore, in order to 
achieve correct physical results, we solve the conformally 
invariant wave equation. We proceed in this manner be- 
cause we are interested in measuring the scalar field far 
from the source and at infinity; in this case it is natural 
to consider future null infinity as the boundary because it 
is the asymptotic boundary of a scalar field propagating 
at the speed of light. 

The paper is organized as follows. In Sec. II we de- 



scribe the construction of boson star space-times. In See. 

III we develop the process of conformal compactification 
for a sphericahy symmetric space-time, that we apply to 
the Schwarzschild and boson star space-times. In Sec. 

IV we present the solution of a massless scalar field on 
these space-times and compare the behavior by solving 
the conformally invariant wave equation. Finally in Sec. 

V we present some conclusions. 



II. BOSON STARS 

Boson stars (BSs) arise from the Lagrangian density of 
a complex scalar field minimally coupled to gravity, that 



For these coordinates the Einstein-Klein-Gordon system 
of equations reads: 
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where kq = 16ttG in units where c = 1, $ is the scalar 
field, the star stands for complex conju gate , and V the 
potential of self-interaction of the field |23, [SO] ■ Notice 
that this Lagrangian density is invariant under the global 
U{1) group, and the associated conserved charge is the 
quantity called the number of particles (defined below) . 
When the action is varied with respect to the metric, Ein- 
stein's equations arise G^^ — kqT^u, where the resulting 
stress-energy tensor reads 
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Boson stars are related to the potential V = tti |$| + 
■| l^l"*. The quantity m is understood as the mass of a bo- 
son and A is the coefficient of a two-body self-interaction 
mean field approximation. The Bianchi identity reduces 
to the Klein-Gordon equation 
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where n$ = 

Boson stars are spherically symmetric solutions to the 
above set of equations under a particular condition: the 
scalar field has harmonic time dependence $(r, i) = 
(/)o(r)e~*"*, where r is the radial spherical coordinate. 
This condition implies that the stress-energy tensor in 
(P|) is time-independent, which implies through Einstein's 
equations that the geometry is also time-independent. 
That is, there is a time-dependent scalar field oscillating 
upon a time-independent geometry whose source is the 
scalar field itself. It is possible to construct solutions for 
boson stars assuming that the metric can be written in 
normal coordinates as 



ds^ 
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The system ([5]) is a set of coupled ordinary differential 
equations to be solved under the conditions of spatial flat- 
ness at the origin a(0) = 1, (^o(O) finite and dr4>o{0) = 
in order to guarantee regularity and spatial fiatness at the 
origin, and (J)q{oo) = in order to ensure asymptotic flat- 
ness at infinity as described in [l^, lU, 123, ImI ; these con- 
ditions reduce the system ^ to an eigenvalue problem for 
a;, that is, for every central value of (pg there is a unique 
uj with which the boundary conditions are satisfied. In 
practice the system of equations is not integrated up to 
infinity, instead the solution is calculated up to a large 
value of r where the boundary conditions are satisfied 
with a given tolerance. Because the system ^ contains 
several constants, it is convenient to rescale the variables 
in such a way that they do not appear. In order to do so, 
the following transformation is convenient: (f>o = \/^(j>Oi 
f = mr, i = Lut, a = —a and A = —2^ Ej. The result 
is that the physical constants vanish from the equations 
and the radial coordinate has units of m and the time has 
units of w. In fact the mass of the boson becomes the 
parameter that fixes the scale of the system. After substi- 
tuting this transformation and removing the tildes from 
everywhere, the resulting Einstein-Klein-Gordon system 
of equations reads: 



a 



dra 
a 



1 



2r 



1 



2" 



ct>l^ + {dr^,f + a'{^l + \Hl) 



r 

9r0O 



^-^-ram{l + \Hl), 



- ci^{l + Hl) 



drCt 

a 

^0 



dra 
a 



a 



= 0. 



(6) 



Notice that the parameter uj now turns into the central 
value of the lapse a due to the rescaling. This is the 
system that is being solved in practice using finite dif- 
ferences with an ordinary integrator (adaptive step-size 
fourth order Runge-Kutta algorithm in the present case) 
and a shooting routine that bisects the value of uj. 
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FIG. 1: Sequences of equilibrium configurations for two val- 
ues of A are shown as a function of the central value of the 
scalar field (/f'o(O); each point in the curves corresponds to a 
solution of the eigenvalue problem and represents a boson star 
configuration. The filled circles indicate the critical solution 
that separates the stable from the unstable branch. Those 
configurations to the left of the maxima represent stable con- 
figurations. The inverted triangles indicate the point at which 
the binding energy is zero. Those configurations between the 
filled circles and the inverted triangles (along each sequence) 
collapse into black holes as a response to a perturbation. Con- 
figurations to the right of the inverted triangles disperse away. 



The solutions of ([6|) define sequences of equilibrium 
configurations like those shown in Fig. [1] Each point 
in the curves corresponds to a boson star solution. In 
each of the curves two important points for each value 
of A are marked: (i) the critical point, marked with a 
filled circle, indicating the threshold between the stable 
and unstable branches of each sequence (that is, config- 
urations to the left of this point are stable and those 
to the right are unstable) as found through the analy- 
sis of perturbations [l9j and full nonlinear evolution of 
the equilibrium solutions |2ll - l23l |: and (ii) the point at 
which the binding energy Eb = M — Nm = marked 
with an inverted filled triangle (see [l^ for this con- 
vention of the binding energy), where N = J j'^cPx ~ 
J §a/— 55'"^[*J^*i9,y$ — ^d„^*](Px is the number of parti- 
cles; that is, the conserved quantity due to the invariance 
under the global C/(l) group of the Lagrangian density 
^. Af = (1 — l/a^)r/2 evaluated at the outermost point 
of the numerical domain is our best guess of the Arnowitt- 
Deser-Misncr (ADM) mass [23; the configurations be- 
tween the instability threshold and the zero binding en- 
ergy point (like configurations 2 and 5) have negative 
binding energy {Eb < 0) and collapse into black holes 
whereas those to the right (like 3 and 6) have positive 
binding energy and disperse away. Those configurations 
to the left of the threshold of instability like 1 and 4, 
that is, stable configurations, possess a negative binding 
energy. Boson stars as black hole mimickers are located 
in the stable branch [9|, llC| , and that is why in this paper 
we only consider stable configurations. 



III. SCRI-FIXING CONFORMAL 
COMPACTIFICATION 

We want to solve for a test scalar field on a fixed 
background geometry corresponding to a boson star and 
compare it with the one calculated on a Schwarzschild 
background. For this we choose hyperboloidal slices that 
show important advantages. At present time, foliations 
with boundaries at future null infinity are being used 
for various app lications, like the solution of perturbation 
equations [32h34| and the study of tails [35| . On the 
other hand, hyperboloidal foliations are useful because 
they reach future null infinity and for asymptotically flat 
space-times, the gravitational radiation and scalar field 
glob al properties are well defined at such a boundary 
j36l |37| . On the other hand, it is possible to compactify 
the space-time in order to work on a numerical domain 
that contains ^7+ , for which it is necessary to compactify 
the spatial coordinates in an appropriate way such that 
the space-time is regular there, at future null infinity. 

Usually, when a spatial coordinate is being compact- 
ified, the metric is singular at the new spatial infinity, 
and the same happens when hyperboloidal slices are com- 
pactified. Among others J38l l39j . a known way of fixing 
this problem is the definition of a conformal metric that 
uses an ad hoc conformal factor that regularizes the sin- 
gular terms resulting from the compactification [23, [S^ . 
The idea behind scri-fixing conformal compactification 
for spherically symmetric space-times is the following: 
(i) use hyperboloidal foliations whose slices reach J'^, 
(ii) compactify the radial coordinate and (iii) regularize 
the metric with an appropriate conformal factor. In or- 
der to proceed, a new time coordinate is introduced as 
t ~ t — h{f), where (i, f) are the original time and space 
coordinates and h{f) is a function called the height func- 
tion [27|, S, 140| . This last transformation has the advan- 
tage that it keeps the time direction invariant. That is, 
regardless of the choice of the h{f), the timelike Killing 
vector has the same form in both coordinate systems. 
A compactifying coordinate is introduced in the form 
f^ — 71, where fl = il(r) is a function that compacti- 
fies r and at the same time is a non-negative conformal 
factor that vanishes at J7+ [23]; in our case we choose 
n^l~rioT the black holes and fl = {l- r'^)/2 for BSs. 
Then J'^ is located at r = 1, where the space-time can 
be rescaled g = fl^g with the conformal factor so that 
the conformal metric g is regular everywhere. 

On the other hand, for the construction of h{f), con- 
stant mean curvature slices are used [40|. The mean ex- 
trinsic curvature k of the initial slice i = is given by 



(7) 



where n'^ is a timclikc unit vector normal to the spatial 
hypcrsurface and positive k means that the slices reach 
future null infinity. In this paper, wc adopt the condition 



that the mean extrinsic curvature k is constant. Such 
condition allows the integration of the equation above. 

We now develop the construction of the conformal com- 
pactification for a spherically symmetric space-time de- 
scribed by the type of metric we start up in the cases 
studied here (Schwarzschild and the boson star): 
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a\f-)dt^ + a\f)df^ + f\de^ + sin^ Odip^), (8) 



where a and a are assumed to be known metric functions, 
and we proceed to construct a hypcrboloidal foliation and 
a scri-fixing compactification. 

By introducing the change of coordinate t = t — h{f) 
the line element takes the form 



d~s^ = -^a^{f)dt^ - 2h'{f)a^{f)dtdf 

+ [a^{f-) - a^{f){h'{f)f] dr^ + f-^{dB^ + sin^ 

where h' = ^. Comparison of this metric with the 3-1-1 
metric ds^ = (-a^ + j'^P^)dt'^ + 2p^^dfdt + j^df^ + 
P{d9'^ + sin^ Odip"^) allows one to read off the gauge and 
metric functions 
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In terms of these functions, the unit normal vector to 
the spatial hypersurfaces pointing to the future is given 
by 
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Given a spacelike slice, we can drag it along the time- 
like Killing vector. This will give a slicing where the time 
translation is along the Killing vector. Now, the mean 
curvature at the initial time slice, given by (O, takes the 
form 
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and can be integrated for constant k: 



k r a{f)a{f)df — C 
where now h!{f) is 
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In this case, it is not easy to find /i in a closed form, so 
in order to have a description of the slices in general one 
has to integrate this function numerically 

On the other hand, in order to perform the scri-fixing 
compactification, we define the compact coordinate r by 
f = j^ and we rescale the original metric by the confor- 
mal factor il. The space-time using scri-fixing conformal 
compactification is finally given by the line element 



ds^ = -a^{r)n^{r)dt^ - 2h' {r)a^ {r){n - rn')dtdr 



+ [a^{r) - a^{r){h'{r)f] (Q - rfl'f 
+r^{de^ +sm'^0dip^), 



dr^ 
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where the functions a{r), a{r) and h'{r) are the functions 
a{r), a{f) and h'{f) evaluated in f = r{r). The conformal 
factor J7 determines various properties of the resulting 
conformal metric. 

The two cases we study here correspond to the 
Schwarzschild and the boson star space-times. In the 
first case, a^ = -^ = {1 — ^^), with which our expres- 
sion ([T4| . using 17 = 1 — r, is reduced to the expression 
for h'{f) obtained by Malcc and Murchadha [401, ^^^ ^he 
final version of the conformally rcscalcd metric reads 
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where 
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The values of k and the integration constant C are re- 
stricted such that P{r) is real. The values of these con- 
stants used in this paper are k = 0.4 and C = 2 in all 
our calculations for the black hole. 

In the Boson Star case: (i) we solve the eigenvalue 
problem (jS]) for a given configuration with a given ADM 
mass; (ii) since such solution is only calculated up to a 
finite -but large- radius, we match the solution to the 
Schwarzschild solution; (iii) then we foliate the resulting 
solution using hypcrboloidal slices with constant mean 
curvature using (jl4[) ; and (iv) finally we construct a con- 
formal metric that contains a compactified spatial coor- 
dinate with the values fc = 3 and C = 0, which we have 
found to be useful for space-times without horizons. 



Since boson star space-times are regular everywhere 
and since we solve the conformally invariant wave equa- 
tion where the Ricci scalar gets involved, it is useful 
to use an adequate confornial factor that allows such 
scalar to be regular at the origin. Following 4l| we use 
O = (1 — r^)/2 which allows a finite value of the Ricci 
scalar at the origin for the Minkowski space-time and we 
found this to be useful in the BS case too. 



TT = Z^dtcf) - IPdrfj) and Tp = 5^0. Then Eq. dH 
transformed into the following system of equations: 






IV. SOLUTION OF THE WAVE EQUATION 

Following our convention, the conformal metric is given 
by ds^ = n^ds^, meaning that the physical metric is the 
one with the tildes. We decide to solve the wave equation 
using the conformal metric because of various reasons: (i) 
the slices are hyperboloidal and reach J'^ at infinity; (ii) 
the spatial coordinate is compactified, and therefore the 
wave function reaches future null infinity at the boundary 
r = 1; (iii) such boundary is null and there is no need to 
impose boundary conditions there. 

In order to use these benefits it is necessary to solve 
the conformally invariant wave equation 



°-6^ 



(j)T{t,r,0,ip) =0, 



(19) 



where R is the Ricci scalar of the conformal metric, D ~ 
V^V;^ also corresponds to the conformal metric and the 
physical scalar field $ is related to the conformal scalar 
field by $ = 0$. In order to study nonradial modes, we 
write separate (t)T{t,r,9,ip) == (f){t,r)Yi,n{0,ip) with Yim 
the spherical harmonics and solve the resulting equations. 
We solve numerically this equation considering a domain 
r G [rmin, !]• For the case of the black holes we choose 
fmin such that it lies inside the event horizon and satisfies 
the need of P{r) being real. For the case of the boson 
stars, since such solutions are real everywhere we choose 



We solve ([191) ^^ ^^ initial value problem using a first 
order variable formulation. In terms of the line element 
([T5)) we can read off again the gauge in terms of the ADM- 
like metric ds'^ = (-d^ + j'^^'^)dP + 2$'y^drdt + j'^dr^ + 
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metric functions 
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which are the final metric functions we use to solve 
the wave equation. We choose to solve p^ as a first- 
order system, for which we define first-order variables 



where we have included the separation of the angular 
part of the scalar field in (|T9| . Finally, we integrate the 
system of equations (|2ip numerically using a finite differ- 
ences approximation on a uniformly discretized domain; 
we use the method of lines with sixth-order accurate sten- 
cils along the spatial direction and a fourth-order Runge- 
Kutta time integrator. 

We use initial data corresponding to an outgoing Gaus- 
sian profile for the scalar field: 



0(0, r) 
V;(0,r) 

7r(0,r) 



-2i!:^0(O,.), 



~^(0,r) 



0(0, r) 



(22) 



At the boundary r = 1 we do not need to impose 
boundary conditions as mentioned above, however at r ~ 
Tmin wc do thc foUowiug: (i) in the case of the black holes, 
since this boundary is located inside the event horizon, 
we do not need to impose boundary conditions given that 
the cones point inward there and we only verify that no 
spurious signals propagate out of the horizon; and (ii) in 
the case of boson stars, we have to deal with the origin 
r = 0, where we first stagger the origin, and second we 



substitute the term A 4^ by 3 
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A. Tests 

We use the solution of the test massless scalar field 
on the Schwarzschild solution in order to have a first 
validation of our implementations. For this we choose the 
case of a Schwarzschild black hole solution with mass and 
show the amplitude of the scalar field in Fig [2j where the 
quasinormal mode and the tail decay for / = 2 are shown 
and validated. From such result we fit the frequency of 
the oscillation which is consistent with the results in 1431 ■ 



B. Solution for Boson Stars 

In order to compare the behavior of the scalar field on 
a black hole space-time with the behavior on a boson star 
space-time we choose boson star configurations that are 
stable, that is, potentially black hole mimickers [l^. We 
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FIG. 2: (Top) We show the absolute value of the scalar field 
amplitude measured at a fixed radius r — lOOOM, for vari- 
ous values of the foliation parameters for a black hole of mass 
M = 1/2, and mode I — 2. We found that the quasinormal 
modes have a frequency uj = 0.749 ± 0.004 which is consis- 
tent with previous studies than indicate to — 0.7473 [42l, and 
the tail decay exponent are both independent of the folia- 
tion parameters. (Bottom) We show the tail decay exponent 
for / = to be of the form --^ t^ measured also at various 
distances from the black hole's horizon. These results are 
consistent with the predictions of measurements by timeliko 
observers (p = — 3 [431) ''■^d null observers at J7^ [p — —2 
[4J]), and also with recent numerical calculations of the poly- 
nomial exponent [2^. The curves represent the measurment 
by detectors located at: infinity, lOOOM, 250Af , llOM, 50M, 
24M and 5M, from top to bottom, respectively. 



also choose two different values of the self-interaction pa- 
rameter A = 0, 20, so that we sample a representative sec- 
tor of the parameter space of boson star configurations. 
For this wc choose the four configurations described in 
Table H 

Wc measure the amplitude of the scalar field at vari- 
ous fixed distances. The absolute value of the scalar field 
measured at various distances is shown in Fig. [3] for the 
four configurations described in the table. It can be no- 
ticed in the figure that for the I = mode, the scalar field 
simply decays and does not show oscillations of any kind, 
pretty much as also found for black holes. However, for 
the mode I = 2 the pair of compact configurations (cases 



TABLE I: Configurations selected to compare the evolution 
of the test scalar field on a boson star space-time. Configu- 
ration A is the most compact one, that is, it has the highest 
2Mgg/rgg ratio, where rgg is the radius of the 2-sphere con- 
taining 99% of the whole mass of the configuration (see [3|); 
in the same sense, configuration B is a very diluted configu- 
ration. The effect of self-interaction is that more mass can be 
maintained by gravity at the price of a more expanded - less 
compact - configuration [31, here configuration C is next to 
the critical configuration and more compact than configura- 
tion D (see Fig. 1). 



A and C) show a set of oscillations, whereas their diluted 
counterparts (cases B and D) do not show a clear train of 
oscillations and the amplitude of the scalar field simply 
drops. This is an important fingerprint, and a property 
to cross-check the boson star parameters of boson mass 
and self-interaction (or compactness) used to construct 
black hole mimickers [3, ll3| • 

Aside from the oscillations, we track the tail-type of 
behavior following the phase of oscillation of the scalar 
field. What has been found for black holes, is - as shown 
in Fig. [5]- that the scalar field amplitude shows a poly- 
nomial tail decay of the type ^ t^ for large values of t, 
where p approaches a constant negative value, and also 
observers located at various positions measure different 
tail decay exponents. We explore the same possibility 
for boson star space-times. Our results for I — arc pre- 
sented in Fig. [H For this we calculate the exponent p 
using the fact that assuming (p{x, t) ~ Kt^ with K a con- 
stant for each position where the amplitude of the field is 

rf log 101 
dlog(t) ' 



measured, the exponent can be expressed by p 
and this exponent may depend on position and time as 
shown for the black hole case in Fig[2| We found in all the 
BS cases that the exponent p does not stabilize around a 
constant value, instead, the exponent remains decreasing 
up to amplitudes of the wave function near to round-off 
error values. For I = 2 we found similar results, although 
the amplitude of the scalar field approaches round-off er- 
ror amplitudes very quickly. 

Finally, we have verified the convergence of our results 
and found that it ranges from fourth to sixth order of con- 
vergence when resolution is increased, as expected from 
the accuracy of our numerical methods. 



V. CONCLUSIONS 

Wc have solved the masslcss test scalar field equation 
on space-times of Schwarzschild black holes and boson 
star configurations. The motivation has been the quest of 
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FIG. 3: (First row) We show the tail decay exponents for two boson star configurations with A = 0. In the first plot we show 
the scalar field amplitude corresponding to the case of central scalar field value 00(0) = 0.25 measured at various distances 
from the center of the star; this boson star configuration is next to the most compact configuration which is the critical solution 
corresponding to 0o(O) — 0.254 [9|; for I = 2 we show a window where oscillations similar to those of quasinormal modes are 
observed, whereas for I — the field simple decays. (Second row) We show the same quantities for the case of a boson star with 
00 (0) = 0.1 which is less compact that the other configuration. (Third and fourth rows) Results corresponding configurations 
with A = 20 and 4>o{0) = 0.1 nearly the most compact for this value of self-interaction and 0o(O) — 0.05 less compact. 
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FIG. 4: We show the polynomial exponent p measured by various observers in time for I = 0. No signs of stabilization of p 
were found. 



restrictions on the boson star parameters acting as black 
hole mimickers. We have found that this test field has a 
very different behavior when evolving on a Schwarzschild 
black hole space-time and when it evolves on a boson star 
background. 

We found oscillations of the test field and signs similar 
to a quasinornial mode ringing for compact configura- 
tions in the I = 2 case, whereas for diluted solutions 
such effect was not as clear. Since we are solving for the 
massless scalar field and not for the perturbation mas- 
ter equation of boson stars - which has to include the 
coupling of the matter the boson star is made of - we 
cannot claim that these are the quasinormal modes of 
boson stars. However we have an indication that not all 
the stable BS configurations show clear oscillations for all 
the potential configurations that may act as black hole 



mimickers. 

Another important sign we have found within the pre- 
cision of our algorithms is that for boson stars, even 
though they correspond to a space-time with curvature, 
no polynomial tail decay was found - at least we did not 
track a stable polynomial exponent. We consider this to 
be an very important fingerprint of gravitational radia- 
tion detected once tails can be measured. 
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